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Q ' Polymer melts represent the flow-history-dependent behavior. To clearly show this behavior, we have 

c/3 , investigated flow behavior of an entangled polymer melt around two cylinders placed in tandem along 

the flow direction in a two dimensional periodic system. In this system, the polymer states around 
a cylinder in downstream side are dififerent from the ones around another cylinder in upstream side 
because the former ones have a memory of a strain experienced when passing around the cylinder 
, in upstream side but the latter ones do not have the memory. Therefore, the shear stress distributions 

P5 ' around two cylinders are found to be different from each other. Moreover, we have found that the av- 

O , eraged flow velocity decreases accordingly with increasing the distance between two cylinders while 

the applied external force is constant. While this behavior is consistent with that of the Newtonian 
fluid, the flow-history-dependent behavior enhances the reduction of the flow resistance. 
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p 

O ■ 1. Introduction 

Predicting flow of entangled polymer melts is difficult because the microscopic states of polymers 
depend on the flow history, namely the history of the experienced strain during the flow. Moreover, 
the quite large number of the degrees of freedom in the entangled polymer melt makes it difficult 
^ ■ to build the constitutive equation that describes the relation among the stress tensor and the strain 

(or strain-rate) tensor in the entangled polymer melt. To overcome the difficulty existing in the flow 
prediction of entangled polymer melts, we have developed a new multiscale simulation [1] composed 
of the macroscopic fluid particle simulation [2, 3] and the microscopic entangled polymer dynamics 
simulation [4-8]. Using the fluid particle simulation where each fluid particle has a polymer simula- 
tor that describes the polymer states in the fluid particle itself, we can accurately consider the flow 
history of the fluid particle and the polymers in the fluid particle. The entangled polymer dynamics 
simulations [4-8] are based on reptation theory [9-11] where the dynamics of a polymer chain is 
constrained in a tube created by the suiTounding polymers because of excluded volume effect and 
entanglements between polymers. Because each polymer chain in the entangled polymer dynamics 
simulation is described with the number of hypothetical entanglement points Z and the tube segments 
r^j{j = 1 , • • • , Z) connecting the entanglement points on the chain, we can decrease considerably large 
amount of the degrees of freedom into a manageable number of the degrees of freedom. The multi- 
scale simulation enables us to simulate the polymer melt flow with remaining the detailed information 
on the entangled states of polymers. 

In Ref. [1], we have considered a flow around a cylindrical obstacle using the new method and 
found that the polymer states are different between the upstream region and the downstream region 
because of the flow-history-dependent behavior of the entangled polymer melt. In that work, we 
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Fig. 1. Schematics of the system. Two infinitely long cylinders are placed in a rectangular system with an 
area L x 2L. The length L - 30ao and the diameter - 6flo of both cylinders are fixed where oo is the 
unit length of the macroscopic fluid simulation. The distance d between two cylinders is a constant parameter. 
The ;c-axis is chosen to be parallel to the flow direction and the y-axis is vertical to the x axis. Because of 
the translational symmetry towards the z-direction perpendicular to the x and y axes, we can focus on the two 
dimmensional flow at the macroscopic level. When d - L, this system corresponds to the one considered in 
Ref. [1]. 



considered the case that there is only one cylinder in the system. If there are two cylinders in the 
system as shown in Fig. 1, the flow history in the upstream region will afi'ect the polymer states in 
the downstream region. Moreover, we expect that the distance d between two cylinders can affect 
the flow behavior of entangled polymer melt. Because the cylinders are placed in tandem along the 
flow direction, the efi'ect of the memory of the flow history can depend on the distance d between the 
cylinders. The purpose of the present paper is to clearly show the flow-history-dependent behavior of 
entangled polymer melt and to investigate the eff'ect of the memory of the flow history using the new 
multiscale simulation [1]. 

2. Multiscale Simulation 

In the multiscale simulation employed here, the polymer melt is described with an amount of 
fluid particles. The position r, and the velocity v,- of the i-th fluid particle are updated according to the 
following equations: 

d7^^'' 

Pi^^-^Pi + ^-o-i+f, (2) 
at 

where / is an external body force with / - {fx,Qi). These equations are integrated with the velocity- 
Verlet algorithm. The spatial derivative of a field variable V/ is calculated with the modified smoothed 
particle hydrodynamics (MSPH) algorithm [2, 3]. The density p, at /-th fluid particle is calculated 
with the usual smoothed particle hydrodynamics (SPH) [12] technique: p, - H jen, "'o^d'"/ - fj\,h) 
where W{r,h) is a Gaussian function with a cutoff length 2h [2]. The pressure pi is a function of 
Pi. The stress tensor cr,- depends on the microscopic states of polymers simulated with the polymer 
dynamics simulation. Using the dual-slip-link model (DSLM) [5, 6] as the microscopic simulation in 
the multiscale simulation, we can obtain the polymeric stress o-p coming from the entangled polymer 
states. DSLM can describe the polymer dynamics larger than the tube diameter b. To include the 
faster dynamics than the tube dynamics (or the dynamics less than the tube diameter), we supply the 
dissipative stress cr^ to the stress tensor o". Because the relaxation time of such faster dynamics 
is negligibly small compared to the tube dynamics, we can regard the dissipative stress as a 
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Newtonian stress tensor with a constant viscosity rj . Then, the stress tensor o",- is assumed to be 
the sum of ctP and a^: 



o-i^o-^ + al (3) 
The polymeric stress (t? is obtained from DSLM as follows: 




where Z is the number of entanglements in a polymer chain, r'. is the y-th tube segment vector on 

the chain, and ctq = {rf lrf^)(TQ is the stress coefficient that transforms the stress unit cr^ in the 
microscopic simulation to that in the macroscopic simulation through the ratio if Irf^"*. The bracket 
(•••), means the average over the polymer chains in the /-th fluid particle. The dissipative stress cr^ is 
assumed to be proportional to the deformation rate tensor D as follows: 

= n^Di, (5) 

where rf is the dissipative viscosity and Di = Vv,- + (Vv;)^- 

The macroscopic fluid particle simulation and the microscopic polymer dynamics simulation 
communicate with each other through the following procedure: 

(1) Ki = (Vv;) is calculated in the fluid particle simulator using MSPH algorithm. 

(2) Ki is transfered to the polymer simulator. 

(3) The polymer states, namely {Z} and {r^}, subjected to /c, are updated for At in the polymer 
simulator. 

(4) cr? is calculated using Eq. 4 

(5) (r? is transfered to the fluid particle simulator. 

The process (3) is composed of several steps, shortly explained as follows (for more information, see 
Ref. [6]): 

(3.1) Subjected to the tube segments {r\} are affinely deformed. 

(3.2) The updated chain length / - T^^^i k' l + si + S2 is relaxed only for a time interval At where 
si and S2 are the lengths of end-segments of the polymer chain out of the tube; The relaxation 
dynamics is composed of the reptation motion and the chain retraction where the equilibrium 
length of the polymer chains, l^q = {Z}eqb, is assumed. 

(3.3) The number of entanglement Z is renewed depending on the length of si (i - 1,2). When 
Si is larger than b, a new tube segment with a length b is created at the end of the chain in 
Si side, and b is subtracted from Sj. At the same time, a new entanglement point is created 
in the other polymer chain randomly selected. This new entanglement point and the new end 
of the polymer tube are coupled with each other, which mimics the entanglement between 
the polymer chains. When .v, is less than 0, the entanglement at the end of the polymer tube 
is erased and the entangled pair is also erased, which mimics the disentanglement process 
(constraint release). 

The procedure from (1) to (5) is the most time-consuming part in the multiscale simulation due to 
the excessively large number of degrees of freedom coming from the large number of fluid particles 
each having the large number of internal degrees of freedom. However, because the procedure from 
(1) to (5) is independent of the other fluid particles, parallel computing drastically improves the time- 
efficiency of the multiscale simulation. 



Table 1. Fixed parameters in this article (cf., Ref. [1]). 



to 




1.0 


Unit time in the macroscopic simulation 


ao 




1.0 


Unit length in the macroscopic simulation 


mo 




1.0 


Unit mass in the macroscopic simulation 


h 




1.0 


Unit time in the microscopic simulation 




i^b<^ ao) 


1.0 


Unit length in the microscopic simulation (Tube diameter) 


0"e 




1.0 


Unit stress in the microscopic simulation 


At 




l.Ox 10-2;o 


Time-integral step 


<Z>eq 




7.0 


Average number of entanglements in equilibrium 






n.5<Tete 


Zero-shear viscosity of entangled polymer melt when <Z)eq = 7 [1] 






200fe 


Disentanglement time (Reptation time) when (Z)eq = 7 [1] 


m 




1.0wo/(aofo) 


Total zero-shear viscosity of the macroscopic flow 


rflm 




0.9 


Viscosity ratio of polymeric viscosity 






0.1 


Viscosity ratio of dissipative viscosity 


fx 




5.0x 10"'*mo/(aofo) 


External body force 


L 




30ao 


System length 


dc 




6.0ao 


Diameter of the cylinders 


Nf 




1,784 


Total number of fluid particles 


Np 




1,000 


Number of polymer chains in a fluid particle 


No 




24 


Number of fluid particles that represent the perimeter of a cylinder 


The cylinder is 


represented with the A^c fluid particles evenly placed on the perimeter. These 



fluid particles fixed on the space also have their own polymer simulator that guarantees the non-shp 
boundary condition assumed on the interface between the polymer melt and the wall [13]. 

The parameters used in the present work are summarized in Table I. Using these parameters, we 
have found that the Reynolds number Re = pUdc/rjo is about 0.3 and the Deborah number De = 
Uzi/dc is about 1.8 where U is the average flow velocity in the system and is set to 0.055ao/?o when 
d = 30ao (cf., Ref. [1]). 

3. Results and Discussions 

Now we investigate entangled polymer melt flows in the system shown in Fig. 1 using the mul- 
tiscale simulation. At first, we focus on two diff'erent cases (a) and (b) where the distances between 
two cyhnders are (a) d - lOao and (b) d - 20ao (cf., Ref [1] when d - 30ao). As is the case in 
Ref. [1], we focus on the magnitudes of the velocity field |v|, the shear deformation rate field \Dxy\, 
and the shear stress field normalized by the total zero-shear viscosity \o-xy\/T]o. The fluid particle data 
are transformed into the data on the square lattice mesh using the hnear interpolation method. Then, 
to decrease the noise of the data, we take an average of the obtained data over the time interval from 
2,OOOfo to 3,OOOfo in the steady state, because the macroscopic states of flow have reached steady 
states before 1, OOO/q in this system. 

The spatial distributions of these data are shown in Fig. 2. Comparing these two cases, as we 
expected, the flow-history-dependent behavior is clearly appeared in the shear stress distribution in 
(a) where the magnitude of the shear stress field around the cylinder in the upstream side is clearly 
difl'erent from that around the other cylinder in the downstream side, while it is difficult to observe 
such a behavior in (b). Moreover, focusing on the velocity fields in (a) and (b), the average flow 
velocity U in (a) seems to be higher than that in (b) because the high velocity region in (a) is broader 
than that in (b). 

To clarify the relationship between U and d, we plot the average flow velocity U against the 
distance d in Fig. 3 (a). For the comparison, the open circles represent the results in the Newtonian 
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Fig. 2. Flow behaviors in cases of the distances (a) d = 10a and (b) d = 20a between two cylinders. The 
magnitudes of the velocity v, the shear deformation rate D^y and the shear stress normalized by the total zero- 
shear viscosity o-„/t]{) are shown here. 



fluid (77*^ = ?7o and if - 0). As shown in Fig. 3, the average flow velocity U decreases with increasing 
the distance d when d < 30ao- Because of the periodicity of the system, the figure shows a symmetry 
on the line d = 30ao, namely the configuration of the cylinders with the distance d for d > 30ao 
corresponds to that for \d - 30ao|. Normalizing U and d a.s U' = {U{d) - U{30ao))/U{30ao) and 
d' = \d - 30ao|/30ao> U' shows a monotonic increasing behavior with d' as shown in Fig. (b). 
The solid line and the dashed line in Fig. (b) are fitting lines obtained using the nonlinear least 
squares fitting. The shapes of these fitting curves are U'{d') = cpd'"^ where cp = 1.411 ± 0.016 and 
ap = 2.979 ± 0.035 for the polymer meh flow and U'{d') = c^^d'"^ where cn = 0.857 ± 0.014 and 
ckn - 2.969 ± 0.051 for the Newtonian flow. Comparing these results, the normalized flow velocity 
U' of the polymer melt is clearly enhanced with increasing d', and the powers of the fitting functions 
for the polymer and Newtonian flows, however, are found to have almost the same value; cp > cn 
and ap - c^N - 3. Namely, the flow-history-dependent behavior enhances the average flow velocity 
(or decreases the viscosity) but does not affect the exponent of the power law function in this flow 
system. 

4. Summary 

We have investigated entangled polymer melt flow around two cylinders in tandem in a two 
dimensional rectangular system, using the new multiscale simulation [1] that is composed of the 
macroscopic fluid particle simulation and the microscopic entangled polymer dynamics simulation. 
In the system, the polymer melt has represented the flow-history-dependent behavior in the shear 
stress distribution. We have found that the flow-history-dependent behavior causes the orientation of 
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Fig. 3. The averaged flow velocity U and the distance d between two cyUnders. The fiUed circles in the graph 
represent the polymer melt flow. For the comparison, the Newtonian flow results (77'' = 770, rf - 0) are also 
shown in the graph using open circles. When d is less than 30a, the averaged flow velocity U decreases with 
increasing the distance d. Because of the periodicity of the system, this graph shows a symmetry on c/ = 30ao. 
namely U(d)foid > 30flo corresponds to J7(60ao-'^)- Normalizing U and (i as U' - (U(d)-U(30ao))/U(30a()) 
and d' - \d - 30aol/30aoj the normalized velocity U' is found to be proportional to d'^ both in the polymer 
melt flow and the Newtonian flow. 



polymer chains, and then reduces the flow resistance. However, the flow history does not affect the 
exponent of the power law function obtained fitting the average flow velocity against the distance 
between cyhnders. 
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